Theoretical evidence for the semi-insulating character of A1N 
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We present ab initio density-functional calculations for acceptors, donors, and native defects in 
aluminum nitride, showing that acceptors are deeper (Be ~ 0.25 eV, Mg ~ 0.45 eV) and less soluble 
than in GaN; at further variance with GaN, both the extrinsic donors SiAi and Cai, and the native 
donor Vn (the anion vacancy) are found to be deep (about 1 to 3 eV below the conduction). We 
thus predict that doped A1N will generally turn out to be semi-insulating in the normally achieved 
Al-rich conditions, in agreement with the known doping difficulties of high-a: Al^Gai-^N alloys. 

PACS numbers : 71.25.Eq, 61.72.Vv, 61.72.Ss 
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III-V nitrides have by now established themselves as a 
key materials system for high-frequency optoelectronics. 
Doping is still a central problem in this area, especially as 
recent developments tend to focus on heterojunction sys- 
tems, and hence on alloys with large In and Al content. 
Doping of either kind has been shown to be increasingly 
difficult as the A1N fraction increases in AL-Gai-aN al- 
loys pL On the other hand, theoretical studies of dop- 
ing have been much less detailed for A1N than for 
GaN |0,p|-||, and practically null for InN. Here we con- 
tribute to the discussion the results of accurate first prin- 
ciple density functional theory calculations for acceptors, 
donors, and vacancies in wurtzite A1N. We find that typ- 
ical candidate acceptors, namely cation-substituting Be 
and Mg, are moderately but sizably deeper than in GaN. 
More interestingly, donors known to be shallow in GaN 
(Si, C, and the nitrogen vacancy) are in fact deep in A1N 
in their substitutional configuration. Our results are in 
general agreement with previous results, when compara- 
ble. An exception is the behavior of donors, which is 
at partial variance with previous studies ||||] , but seems 
to agree with experimental data indicating serious n- 
type doping difficulties for A1N and high- a; Al^Ga^^N 
alloys. 

Method We use local-density-functional-theory 
(LDA) H ultrasoft-pseudopotential plane-wave calcula- 
tions of total energies and forces in doped A1N wurtzite 
supercells typically encompassing 32 atoms, and with 
plane wave cutoff of 25 Ry (for further technical de- 
tails see Refs. ||[7 10 ), to predict from first principles 



the formation energies and thermal ionization energies 
of cation-substituting Be, Mg, Si, and C, and of the N 
and Al vacancies in A1N. The carriers concentration at 
temperature T due to an impurity or defect with thermal 
ionization energy e, and a formation energy Ef olm , is 



N s exp (--Eform/kBTg) exp(-e/k B T), 



(1) 



asssuming impurity incorporation or defect creation in 
thermal equilibrium at a growth temperature of T g , and 
with A s =2.44xl0 22 cm -3 available cation or anion sites 



(i.e. half the theoretical atomic density of A1N). Thus, 
the largest the formation and ionization energies, the less 
efficient the doping. A non-zero formation entropy (ne- 
glected here) will of course enhance the dopant concen- 
tration. The formation energy for an impurity in charge 
state Q is 



-Eform(Q) — E 



X 



+ Q(jjL e + E$), (2) 



with fi e the electron chemical potential (synonimous with 
the Fermi energy Ep in our T=0 calculations), E tot (Q) 
the total energy of the fully-relaxed defected supercell in 
charge state Q, E® its top valence band energy, n x and 
/x x the number of atoms of the involved species (X=A1, 
N, impurity) and their chemical potentials. The latter 
potentials are determined by the equilibrium conditions 
with A1N and the compounds (if any) of the specific 
impurity with Al or N. The structures and formation 
enthalpies of the solubility limiting compounds (AI2O3, 
SiaNi, Be3N2, Mg3N2, and d— C) are calculated ab ini- 
tio. We generally assume the highest ^ lm P urlt y compati- 
ble with the relevant solubility limit. Then, only only one 
independent chemical potential is left (e.g. /i N ), and its 
value is determined by the imposed (N-rich <-> Al-rich) 
growth conditions. 

The thermal ionization energy e[0/— ] of a single ac- 
ceptor is by definition the formation-energy difference of 
the two charge states Q=0 and Q=-l at /i e =0, and it 
corresponds to hole release from (i.e., electron promo- 
tion into) the acceptor extrinsic state with concurrent 
geometric relaxation. For a single donor, the analogous 
quantity is e[+/0], the formation-energy difference of the 
charge states Q=+l and Q=0 at /i e =0. This quantity is 
the distance of the energy level from the valence band. 
The thermal ionization energy of the donor electron into 
the conduction band is E gap — e[+/0]. Thus, the theoret- 
ical evaluation of donor ionization energies with respect 
to the conduction edge requires a decent estimate of the 
fundamental gap. The DFT-LDA eigenvalue gap is in- 
appropriate for two reasons: first, it is notoriously in- 
accurate in general fli]]; second, it seems inconsistent to 
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compare the DFT-LDA eigenvalue-difference gap at fixed 
number of electrons, with our ionization levels extracted 
from total energy differences of variable-electron-number 
systems. 

We thus proceed to evaluate the gap in analogy to im- 
purity levels, specifically as the ionization energy e[Q/— ] 
of a defectless cell, i.e. the difference in total energy be- 
tween bulk supercells with N and N+l electrons respec- 
tively. This way of proceeding is technically consistent 
with the treatment of the defects, and the final result is 
affected by essentially the same systematic errors: our 
result of 5.6 eV for the gap, within about 10 % of the 
experimental value of 6.2 eV, indicates that the correct 
theoretical gap to be compared to donor levels is indeed 
somewhere near the experimental value. Yet another, 
conceptually more satisfactory estimate would be that 
dictated by ASCF theory (see e.g. Ref. @), whereby 
the gap is the difference e[0/— ]- e[+/0] for the defectless 
cell. Using this expression, we get 6.0 eV, even closer to 
experiment. The gap calculated as discussed above, and 
the ensuing depth of the impurity levels, should be taken 
as qualitative indications ]l3[ . 

Results - The results of our calculation are summa- 
rized in Figure |l|, for both Al-rich (upper panel) and N- 
rich (lower panel) growth conditions. The thermal levels 
of acceptors, referred to the valence band top, are 0.25 
eV for BeAi and 0.45 eV for MgAi- Thus, these candi- 
date acceptors are both appreciably deeper than in GaN, 
where the levels of BcQa and MgQ a are at about 0.1 and 
0.2 eV, respectively 
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FIG. 1. Formation energy of impurities and native defects 
in A1N, for N-rich and Al-rich growth conditions. Arrows in- 
dicate the best Fermi levels achievable in typical conditions 
(T 9 ~ 1200 K). 

The extrinsic cation-substitutional donors are found to 
be deep, with thermal levels at 0.9 eV (Si) and 1.2 eV (C) 
below the conduction band (the position of the latter is 
the valence band top position plus the gap estimated as 
discussed previously). Si and C are stable in the substi- 
tutional configuration, and are found to be deep therein 
JLlj . Inspection of the partial charge density of the im- 
purity state of SiAi confirms its strong localization. For 
both C and Si, the neutral state is never stable, i.e. both 
are negative-ZJ centers. The calculated U's are -0.22 eV 
and -0.17 eV for Si and C respectively. 

As for the vacancies, Vai is a multiple acceptor with 
a formation energy diving to zero and below in n-type 
conditions. This relatively unsurprising results parallels 
closely that for the Ga vacancy in GaN || , and the pre- 
vious results on A1N by Mattila and Nieminen Q. The 
nitrogen vacancy is a donor in p-type conditions, as in 
GaN: however, its first thermal level is very deep, indeed 
close to midgap. Negative states of Vn are realized up 
to -3, in agreement with previous results Q. We did not 
investigate the +2 and +3 states of Vn found in Ref. ||, 
but their existance does not affect our conclusions. The 
realization of high charge states is not uncommon in in- 
sulators and wide gap semiconductors, a known example 
being oxygen vacancies in a— quartz p5| . 

The implications of these results for doping are eas- 
ily analyzed using the theoretical framework by van de 
Walle et al. fll6| |. For our present purposes we just need 
to recall that in the final analysys, in the presence of 
competing donor and acceptor species, the Fermi level 
turns out to be pinned at about the intersection of the 
lowest donor and acceptor formation energy curves. The 
Fermi levels obtained respectively via n and p doping for 
typical realistic conditions (T g ~ 1200 K, T oper ~ 300 
K) are denoted by arrows in Figure |l|. Doping is clearly 
very inefficient in all cases; let us examine the reasons for 
that. 

As for p-doping, it is clear that in Al-rich conditions 
the N vacancy compensates the acceptors, pinning the 
Fermi level at about 2-2.5 eV above the valence band. 
In the more favorable N-rich conditions, this compensa- 
tion is not effective (even if the 3+ state of Ref. || is 
accounted for). However, the formation and ionization 
energies of the acceptors are sizably larger than in GaN, 
and both the dopant and the carrier concentration will 
be substantially lower than in GaN, making p-doping of 
A1N even less efficient than in GaN. The shallowest ac- 
ceptor will produce in the most favorable case a Fermi 
level at ~ 0.5 — 0.6 eV above the valence (BeAi in N- 
rich conditions at a MOCVD-like growth temperature 
T g ~ 1000 — 1200 K). The co-incorporation mechanisms 
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involving H [?]]|] or O (?|,[l7j found to be effective in GaN 
remain to be investigated in A1N. 

As for n-doping the situation is even worse. Firstly, 
the shallowest level found here, that of Siai (the natural 
choice, as Sica is the standard shallow donor in GaN), 
is about 1 eV below the conduction band (C has a huge 
formation energy and is always unfavorable). Secondly, 
vacancy compensation is always in action. The nitrogen 
vacancy will compensate Si in Al-rich conditions, and the 
Al vacancy will in N-rich conditions. Thus, in no event 
will the Fermi level be higher than about 2 eV below 
the conduction band, in both N-rich and Al-rich growth 
conditions, n-doping of A1N will therefore be especially 
difficult, if at all possible. 

In summary, ab initio calculations predict that accep- 
tors and donors in use for GaN are not suitable for A1N 
and high- a; A^Gai^N alloys. The acceptors Be and Mg 
have high formation energies and fairly high ionization 
energies, and suffer from nitrogen vacancy compensation 
in Al-rich conditions. The donors Si and C are deeper 
than ~ 1 eV below the conduction band, and therefore 
unusable for doping of A1N. In addition, in N- and Al-rich 
growth conditions, respectively, the Al vacancy and the 
N vacancy will act to compensate the donors. While the 
acceptor behavior of the Al vacancy is unsurprising, the 
deep donor character of the N vacancy is at variance with 
the results for GaN. In this respect our findings support 
the previous results by Mattila and Nieminen [Q over 
those of Stampfl and van de Walle j|. While moderate 
p-type doping might be obtained in N-rich conditions, in 
Al-rich conditions A1N will be always semiinsulating. 
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